Competition between Intra-community and Inter-community Synchronization 
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In this paper the effects of external links on the synchronization performance of community 
networks, especially on the competition between individual community and the whole network, are 
studied in detail. The study is organized from two aspects: the number or portion of external 
links and the connecting strategy of external links between different communities. It is found that 
increasing the number of external links will enhance the global synchronizability but degrade the 
synchronization performance of individual community before some critical point. After that the 
individual community will synchronize better and better as part of the whole network because the 
community structure is not so prominent. Among various connection strategies, connecting nodes 
belonging to different communities randomly rather than connecting nodes with larger degrees is 
the most efficient way to enhance global synchronization of the network. However, a preferential 
connection scheme linking most of the hubs from the communities will allow rather efficient global 
synchronization while maintaining strong dynamical clustering of the communities. Interestingly, 
the observations are found to be relevant in a realistic network of cat cortex. The synchronization 
state is just at the critical point, which shows a reasonable combination of segregated function in 
individual communities and coordination among them. Our work sheds light on principles underlying 
the emergence of modular architectures in real network systems and provides guidance for the 
manipulation of synchronization in community networks. 

PACS numbers: 05.45.Xt, 87.18. Sn, 89.75.Fb 



I. INTRODUCTION 



Communities exist ubiquitously in all kinds of net- 
works [li^tSl and synchronization on community networks 
is of great importance in the biological and social net- 
works. It has been shown that the sparsity of connec- 
tions between different communities hinders the global 
synchronization of complex networks 0-01 • Factors that 
affect the synchronization of community networks have 
been intensively studied and strategies that can achieve 
global synchronization have been proposed [tUToI]. Syn- 
chronization transition process of community networks 
have also been studied [ll|, [l^ . Utilizing synchronization 
to detect community structure has attracted a great deal 
of attention recently p"3 - [l7j . Moreover, the dynamical 
modules of networks with and without clear communities 
have also been studied in detail [l8l - (23 |. Very recently, 
a class of small- world networks with spatial and network 
modularity were obtained by evolving the arrangement of 
nodes in space and their corresponding network topology 
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[23|. Besides, some other interesting topics, such as the 
synchronization interfaces and overlapping communities 
in complex networks, were also studied |24| . 

In the study of synchronization on community net- 
works, how the competition and coordinations of the 
links within community and between communities shape 
the synchronization properties are still not clearly under- 
stood. In this paper, we focus on the effects of external 
links (links connecting nodes belong to different commu- 
nities) on the synchronization performance of commu- 
nity networks, especially on the individual community. 
We study (i) how the synchronization performances of 
global network and individual community are affected by 
the changing of the number of external links, (ii) whether 
connecting nodes with larger degrees of different commu- 
nities will ensure better global network synchronizability, 
and (iii) the synchronization property of the cortical net- 
work of cat in terms of the above two aspects. This work 
will shed light on the understanding of the impacts of 
the number and the connecting strategy of external links 
on the global network synchronizability and individual 
community synchronization performance. This work is 
organized as follows: in the next section, we introduce 
the dynamical equation of each node and the community 
network model. We will present the two aspects of stud- 
ies on the community network models in the sections III 
and IV, and the study on cat cortex in section V. Then 
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we will give our discussion and conclusion in the last sec- 
tion. 



II. DYNAMICAL EQUATION AND 
COMMUNITY NETWORK MODEL 

To investigate the synchronization on complex net- 
works, dynamical systems are often taken as nodes and 
the couplings between different systems are the links, 
thus the dynamical equation of each node in a complex 
network with N nodes is 

N 

i.=F(x,)-^5]G,,H(x,), i^l,...,N, (1) 

where x = F(x) is the individual dynamics, a is the 
overall coupling strength, (fc) is the average degree of the 
network, H(x) is the output function and Gij is the el- 
ement of coupling matrix. In our simulations, if node i 
is coupled by node j, Gij < 0, otherwise Gij — 0, and 

the diagonal element Ga = — J2f=i j^i ^ij to ensure the 
sum of the elements in a row is so as to make the syn- 
chronization manifold an invariant manifold. In many 
previous work, the global network synchronizability were 
measured by the eigenvalues of the coupling matrix with- 
out calculating the iteration of oscillators on the networks 
according to the master stability function [2^ [2^ . In this 
paper, we study not only the global network synchroniz- 
ability but also the individual community synchroniza- 
tion performance, thus the master stability function is 
insufficient to our study, so we take Rossler oscillator as 
dynamical node to fulfill our simulation: 

X -y-z, 

y = x + ay, (2) 

i = 6 + z{x — c), 

where a = 0.2, h = 0.2 and c ~ 7.0. We take the output 
function H(x) ~ (a;, 0,0). In section II and III, all the 
links are un-weighted and the coupling matrix is Lapla- 
cian matrix, i.e., Gij = — 1 (node i and j are connected) 
or Gij = (node i and j are not connected) for off- 
diagonal elements, and Ga — ki for diagonal elements. 
The correlation between each pair of nodes Cy is calcu- 
lated to measure the synchronization performance, and 
the average correlation of all pairs of nodes in the whole 
network Gn and that of pairs in a community Cc are 
taken as the measure of the synchronization performance 
of global network and individual community, respectively. 

To implement our study effectively, we design the fol- 
lowing community network model: take m networks as 
subnetworks to generate a community network and these 
subnetworks will be the communities, then rewire the 
internal links (links connecting nodes within the same 
community) to external links with some strategies, which 
will be discussed in detail. When more and more exter- 
nal links emerge, the modularity will decrease gradually 
till the network becomes a homogeneous one. 
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FIG. 1: (Color on line) The synchronization performance of 
global network Cat (a), of individual community Cc (b), and 
dynamical modularity Qd (c) as functions of the coupling 
strength a at different topological modularity Qt (or external 
link number Le)- (d) Cjv, Cc and Qd at very small coupling 
strength for Qt ~ 0.500. There are 2 communities in the 
network, each one is a BA scale-free network with 100 nodes, 
and the average degree in each community is 16. Each plot is 
obtained after the averaging over 50 network configurations 
and 10 initial states of each configuration. 

We follow Ref. [13, [HI to quantify modularity of a 
network with m communities as 

L 

Q^J2(eu-af), (3) 
I 

where a; — eik and eik is the fraction of to- 
tal strength of the links in the network that connect 
the nodes between the communities I and fc, namely 
eik = i^/W)J2iei,jek^'.j^ '^ith being the connec- 
tion strength between two nodes and W the total of Wij 
in the network. The topological modularity Qt measure 
the strength of community structure of community net- 
works. Synchronization clustering can be quantified by 
the dynamical modularity Qd computed by taking Wij 
as the dynamical similarity between two nodes. In par- 
ticular, Wij = (1 -|- Cij)/2.0, where Cij is the correlation 
between pair of nodes i and j. 

III. EFFECTS OF EXTERNAL LINK NUMBER 

To study the effects of the number of external links Le 
on the synchronization of global network and individ- 
ual community, we take several BA scale-free networks 
as subnetworks, and construct community networks with 
the following operations: (i) random select an internal 
link, (ii) cut one end of it and (iii) rewire it to a random 
selected node in the other subnetworks, then repeat the 
operations (i) to (iii) till the desired community network 
is obtained. This connecting strategy is called random 
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FIG. 2: (Color on line) The change of state x with time t 
at couphng strength (a) 0.014, (b) 0.026 and (c) 0.060. Two 
black curves and two red curves present nodes belong to the 
two different communities separately and the two solid curves 
present nodes connected by external link and the two dash 
curves present two nodes which are neither connected to each 
other nor to the other two nodes, (d) The change of d_B 
(square line) and di (circle line) with the coupling strength. 
The topological modularity of the adopted network is 0.500 
and the other network parameters are the same as in Fig. [1] 



connection. We have investigated networks composed by 
different number of communities, and found no essen- 
tial differences on the synchronization property. Thus, 
in our simulation we take the simplest community net- 
work model with two communities. 

In Figs. [l{a) and (b), we presented the changes of Cn 
and Cc with coupling strength a at different topological 
modularity Qt (corresponding to different number of ex- 
ternal links). Clearly, at the same coupling strength, with 
the decreasing of topological modularity, the global net- 
work synchronizability will be better and better, which 
is consistent with previous studies jj^iGj. However, the 
change of global synchronization performance respects to 
the coupling strength is of great difference: all the curves 
first increase sharply at very small coupling strength [re- 
gion I in Fig. [TJa)] and then increase slowly in a broad re- 
gion of coupling strength [region II] . With the further in- 
creasing of coupling strength, networks with fewer exter- 
nal links (larger Qt) will show degraded synchronization 
performance, while those with more external links will 
be much better [region III] . After that all the curves will 
reach steady states [region IV] . As for the synchronization 
of individual community, Cc will increase monotonically 
with the increasing of coupling strength. However, Cc 
does not show a monotonic dependence on the number of 
external links, as seen in Fig. [Ijb). We will investigate 
the change of Cc with respect to Qt in more detail later 
in this paper. We also study the effects of external links 
on dynamical modularity Q d , which is obtained from the 
correlation matrix (cij). Figure [TJc) shows that dynam- 
ical modularity Qd has two maxima for each curve as 
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FIG. 3; (Color on hue) The change of Cn (a), Cc (b) and 
Qd (c) as with respect to Qt at various coupling strength a 
as indicated. The other network parameters are the same as 
in Fig. [U 



a function of the coupling strength, one in region I at a 
very small coupling strength and the other in region HI. 

In the following we investigate into the mechanism of 
these phenomena. As seen in Fig. [Ud) for the region 
of small couplings, Cc grows faster than Cat, leading 
to enhanced dynamical modularity. Then Cn becomes 
large too and the dynamical clusters become smeared 
and is reduced. It worth mentioning that very small 
coupling strength could make the network go into very 
good phase synchronization state, and at the coupling 
strength a = 0.030 the initial isolated subnetworks will 
reach complete synchronization state. To show explicitly 
the synchronization state of nodes in different communi- 
ties, we plotted the time series of the variable x{t). In 
each community two nodes (not directly connected) are 
selected, one having some external links and the other 
one only having internal neighbors. By comparing the 
four curves we found that phase synchronization through- 
out the networks occurs at very small coupling strength 
even for very few external links, as seen in Fig. HJa). This 
corresponds to the slowly increasing region of Cn [region 
II in Fig. 1(a)]. With the further increase of the coupling 
strength, the networks move to a regime of fierce com- 
petition: for the networks with very few external links, 
individual community get ready to reach complete syn- 
chronization state and the synchronization of nodes in 
the same community is rather good, but because of the 
random initial states different communities oscillate at 
different phase, and the small number of external links 
fail to bring the phases of the two communities to each 
other, thus the global synchronization shows even worse 
performance compared to weaker coupling strength. This 
phenomenon can be seen in Fig. [2][b) . The good synchro- 
nization performance within individual communities and 
weak synchronization between them result in the promi- 
nent dynamical modularity (region HI, Fig. l(a, c)). As 
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for the networks with enough external hnks, communi- 
ties yield to the effects of external links, and there is 
no significant difference between intra-community syn- 
chronization and inter-community synchronization [Fig. 
[2jc)] , leading to reduced dynamical modularity Q d (Fig- 
1(c)). And when the number of external links is between 
these two situations, the individual community tries to 
oscillate independently but is inevitably affected by the 
others to some extent, which results in degraded syn- 
chronization performance of individual community. This 
is the reason for the decreasing of Cc at some topological 
modularity Qt (Fig. 1(b)). When the coupling strength 
becomes stronger and stronger, the state information ex- 
changes more smoothly between different communities 
and the states of different communities oscillate almost 
fully synchronously (Figure (^fc) ) . The global synchro- 
nization increases again and the dynamical communities 
become not so prominent. This procedure of competition 
and coordination can be further manifested by the state 
differences of a node to its neighbors in the same com- 
munity [di) and in the other communities (de). Here 
the difference for a node i and its neighbors j in a com- 
munity S is defined as di = ^ {{xj^s — Xi^), where xj^s 
is the average state of I's neighbors in community S and 
(•) denotes averaging over time, dj and ds averaged over 
nodes are shown as functions of connection strength a 
in Fig. md). It is evident that dj almost keeps decreas- 
ing with cr, which indicates that the state of a node gets 
close to its neighbors in the same community gradually. 
However, the curve of ds has a distinct peak at about 
a* = 0.026, which shows that in this region of coupling 
strength the state of a node is far away from its neighbor 
in the other community. In summary, we can see that at 
very weak coupling strength, the external links could ef- 
fectively transmit the state information of different com- 
munities to make the nodes of the whole network have 
the similar phase, and at very strong coupling strength, 
the external links could also be efficient to bring together 
the states of different communities. However, when the 
coupling strength is between these two cases, the syn- 
chronization in communities are too strong individually 
and the external links fail to synchronize them. 

We also investigate Cn, Cc and Qd as functions of the 
topological modularity Qt with fixed coupling strengths 
(Fig. [3]). Obviously, the synchronization performance of 
global networks is degraded with the increase of Qt, and 
dynamical modularity Qd will keep about zero at small 
Qt and increase sharply when the community structure 
becomes prominent enough at large Qt- Interesting be- 
havior happens for synchronization of individual com- 
munity where Cc display a pronounced minimum. The 
change of Cc can be understood as follows: when inter- 
nal links are rewired to external links in the beginning 
[Qt decreasing), the number of internal links is reduced 
and the number of external links is increased, synchro- 
nization within community becomes weaker. In fact, only 
decreasing the number of internal links or only increas- 
ing the number of external links will both reduce Cc- 
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FIG. 4: (Color on line) The change of Cc, Cn and Qd with 
a at Le =64 and 196 for different rewiring strategies. There 
are 2 communities in the network, each one is a BA network 
with 100 nodes, and the average degree in each community is 
16. Each plot is obtained after the averaging of 50 network 
configurations and 10 initial states of each configuration. 



At this stage, the individual communities still oscillate 
largely independently {Cc ^ 1 > Cn)- When there are 
enough external links, as part of the whole network the 
individual communities are not prominent at all and it is 
no longer independent and Cc ~ Cn, corresponding to 
Qd-0- 

Therefore, the minimum of Cc corresponds to a criti- 
cal point Qtc- When Qt > Qtc, communities oscillate 
rather independently and the communication between 
communities are weak, and for Qt < Qtc , the communi- 
cation between communities are smooth but the dynam- 
ical independency of individual community is not promi- 
nent. So the critical point represent a balance region 
where the individual community is still clearly indepen- 
dent while the information transmission remains effective 
between different communities. 



IV. EFFECTS OF CONNECTING STRATEGY 
OF EXTERNAL LINKS 

Intuitively, not only the number of external links but 
also the connecting strategy of external links could affect 
the synchronization performance of community networks. 
For example, connecting nodes from different communi- 
ties with larger degrees seems to be more efficient for 
information exchange between the communities, there- 
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fore, might be superior for global synchronization. In 
this section we design several connecting strategies to in- 
vestigate whether it will be better for global network syn- 
chronization to connect nodes with larger degree rather 
than nodes with smaller degree. 

To achieve our aim, we also take two BA scale-free 
networks as communities and take the following connect- 
ing strategies: (i) keep one end of external links on some 
nodes with largest degree (big nodes) and connect the 
nodes with randomly selected nodes in other communi- 
ties (hub-random connection), (ii) connect nodes selected 
in different communities with probability p oc kf, where 
ki is the degree of node belongs to the community (prefer- 
ential connection), and (iii) connect big nodes belong to 
different communities (hub-hub connection). Compared 
to the random connection, these three connecting strate- 
gies make the distribution of external links biased more 
and more to some big nodes. The selection probabil- 
ity in preferential connection strategy guarantees that 
the big nodes of different communities will not only con- 
nect densely among themselves, but also connect to some 
smaller nodes. Such a connection scheme is consistent 
with the structure of brain cortex network that the hub 
areas are almost fully connected among themselves to 
form a hyper-community overlapping different communi- 
ties [1^. Moreover, to compare the effects of degree of 
connected nodes, we also designed the middle/ small de- 
gree connection strategy, which is characterized by con- 
necting nodes with middle / small degree in different com- 
munities. For all the strategies, once the external link 
number is fixed, the topological modularity is almost 
the same, but the synchronization performances for net- 
works connected with different strategies are obviously 
different. Figure 3] shows the comparison results of these 
strategies. From the figure we can conclude that the dis- 
tribution of external links homogeneously on nodes with 
different degree will lead to the best global synchroniza- 
tion performance but the worst dynamical modularity. 
If the external links must concentrate on several hubs, 
connecting large nodes will be much better. 
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FIG. 5: (Color online) (a) The cortical network of cat. A 
node represents a functional region of the cortex and a link 
represents the existence of fiber projection between two re- 
gions. The different symbols represent different connection 
weight: 1 (o sparse), 2 (• intermediate) and 3 (* dense). The 
communities/functional subdivision V, A, SM and FL of the 
network are indicated by the solid lines, (b) The change of 
synchronization performance of global network (Cjv) and of 
the four communities (Cv, Ca, Csm and Cfl) vs. the cou- 
pling strength a. Each plot is obtained after the averaging of 
10 initial states. 



nization properties of the preferential connection scheme 
as in the real cat cortical network allow quite strong 
global synchronization and meanwhile maintain interme- 
diate dynamical modularity, again pointing to a mean- 
ingful balance between segregation and integration. 

In order to study the relevance of this interesting ob- 
servation, we study dynamics of realistic cortical neural 
network of cat. Cat cortical network [37]] [Fig. [5] (a)] has 
53 nods, 826 weighted links and four communities that 
carry out the four functions: visual (V), auditory (A), 
somato- motor (SM) and frontolimbic (FL). In this net- 
work, each node represents a brain area which are com- 
posed of huge number of interacting neurons. The rhyth- 
mic activity of such neural ensemble can be modelle d by 
neurophysiologically realistic neural mass oscillators [3g] . 
The dynamical equation of coupled neural mass model is 

M 



V. SYNCHRONIZATION IN MODULAR 
CORTICAL NETWORKS 

Modular organization is prominent and synchroniza- 
tion dynamics is of special important for functioning in 
neural systems [30l - [34| . 

The critical competition regime in the synchronization 
of community network is of special importance from the 
viewpoint of information processing, where the special- 
ized processing using the dynamical communities can co- 
exist with the global exchange of information and inte- 
gration from specialized communities. In neural systems 
a combination of these two ingredients, the functional 
segregation and integration, is believed to be crucial to 
underly the structural and dynamical organization for 
effcient and diverse functioning (ssl |36| . The synchro- 



Aaf(vf -vl)-2avf -a^vf, 
BbC4f(C3v[) ~ 2hv\ - h'vl. 



(4) 



Aa 



N 



C2f(CivP)+p,(t) 



(S) 



Y.W^,f(v 



- 2avf-a^vf, 



where vf , vj and vf are the average post-synaptic 
membrane potentials of pyramidal cells, excitatory in- 
terneurons and inhibitory interneurons of the area i, re- 
spectively. A static nonlinear sigmoid function f(v) = 
2eo/(l + e'^(''°-^)) converts the average membrane poten- 
tial V into an average pulse density of action potentials. 
Here Vq is the post-synaptic potential corresponding to a 
firing rate of Cq, and r is the steepness of the activation. 
The parameters A and B represent the average synaptic 



6 



0=20 



0.29 


V 






• 


(a) 


0.24 


0.28 


\ 

• 








\ 




0.22 


0.27 


« 
• 


• 






\ 

• 


• 


0.20 


0.26 


• 


* 










0.18 



SM 



FL 



(b) 0.32 
0.31 
0.30 



•. 0.29 



(c) 



if 



0.26 ■ 

0.25 

0.24 



(d) 



0.0 0.2 0.4 



0=30 

O 0.38 



I\ f 

I M 

°H.O 0.2 0.4 



0.47 

=40 0-46 

O 



0.45 
0.44 



0.0 0.2 0.4 



0.0 0.2 0.4 

0.32 /■•-.^ (f) 

0.30 *^ 

0.28 

0.0 0.2 0.4 
0.36 



0.0 0.2 0.4 0.0 0.2 0.4 
n0.40| — 



0.40 
0.38 
0.36 



(g) 



(h) 



(j) 0.50 



0.34 
0.32 



/ • ^]0.48 



46 
0.44 



0.2 0.4 


0.3g 


0.2 0.4 




0.46 














*. f 




0.44 
















' \i 


0.42 





0.0 0.2 0.4 



0.0 0.2 0.4 



0.0 0.2 0.4 



FIG. 6: (Color on line) The change of Cc vs. Qt at a =20.0, 
30.0 and 40.0. The stars present the results of the original 
cat cortical network. Each other plot is obtained after the 
averaging over 50 network configurations and 10 initial states 
of each configuration. 



gains, 1/a and 1/6 the average dendritic-membrane time 
constants. Ci and C2, C3 and C4 are the average num- 
ber of synaptic contacts, for the excitatory and inhibitory 
synapses, respectively, pi represents afferent inputs from 
subcortical systems. A more detailed interpretation and 
the standard parameter values of this model can be found 
in [38i|. We take the parameters as in Q so that the 
model generates alpha band periodic oscillations. As in 
reference [s^, in our simulations we take the subcortical 
input as Pi(t) = PQ + ^i(t) where S,i{t) is a Gaussian white 
noise with standard deviation D = 2. Wij is the coupling 
strength from area j to area i. We normalize coupling 
strength cr by the mean intensity (S) where the connec- 
tion intensity Si = the total input weight to 
node i. 

In this system of noisy oscillations, we also use the 
correlation between nodes to measure the network syn- 
chronization performance. Fig. [5] (b) shows the change 
of Cat and C'c of the four communities with the coupling 
strength a. Unlike in Fig. 1(a), here we do not observe 
the region III where Cn decreases with a; however, there 
is a regime at intermediate couplings where C n increases 
slowly. 

In order to see if the real system could work in the 
balance/critical states, we compare the synchronization 
performance of the real network with rewired networks 
having larger or smaller topological modularity. This is 
obtained by rewiring the internal links to external links or 



rewiring the external links to internal links of the original 
cat cortical network. Figure [S] shows the change of Cc of 
the four communities with respect to Qt- Interestingly, 
in a broad coupling strength region, the original network 
(star) is just at the critical point of topological modu- 
larity where Cc is a minimum simultaneously in several 
or all the communities. These results confirmed our ex- 
pectation that the real cortical network is such organized 
that it allows efficient global integration under the con- 
dition that the functional segregation is also maintained 
simultaneously. This finding provides a detailed mecha- 
nism for our previous observation [22| that the dynamical 
complexity measuring a balance between segregation and 
integration is optimal in cat cortical networks. 



VI. DISCUSSION AND CONCLUSION 

In summary, we have studied in detail the effects of ex- 
ternal links on the synchronization performance of com- 
munity networks. We have revealed an interesting com- 
petition between synchronization of the global network 
and the individual communities. With the increasing of 
the number of external links the global network synchro- 
nization will be enhanced but the synchronization per- 
formance of individual community will be degraded till a 
critical point where the community structure is no longer 
prominent. Afterwards, synchronization within the com- 
munity increases again as part of the global network. We 
also investigated the impact of various connection strate- 
gies on the global and community synchronization. We 
showed that connecting nodes selected randomly in dif- 
ferent communities will ensure better global network syn- 
chronization, but weak dynamical modularity. The dis- 
tribution of the external links mainly among the hubs 
nodes would allow both strong global synchronization 
and clear dynamical modularity. 

Interestingly, these discoveries in generic models are 
demonstrated to be relevant in a realistic cat cortical net- 
work with simulated neural population activities. Com- 
paring the synchronization properties to rewired net- 
works with larger or small modularity, we found that the 
real network is just at the critical point. The hubs in this 
network are also responsible for inter-community connec- 
tions similar to the preferential connection schemes in the 
model. Our analysis indicates that the cat cortical net- 
work is such organized that it allows both segregated per- 
formance with the communities and efficient integration 
of the whole network. 

Our work has provided a deeper understanding how 
the external link number and connecting strategy affect 
the synchronization of community networks as a whole 
as well as individual community. These results not only 
present the possible reason for the real cortical network 
to evolve to the current community structure from the 
dynamical point of view but also provide useful methods 
to regulate the synchronization of community networks 
for potential applications. 
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